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Abstract 



Results for the energy dependence of the elementary kaon production cross sections 
q\ ' in proton-proton collisions are reported. Calculations are performed within an extended 

version of the resonance model which was used for the previous studies of elementary 
kaon production in pion-nucleon and pion-A collisions. Although the model treatment is 
within the empirical tree level (observed widths for the resonances are used), it is fully 
relativistic, and includes all relevant baryon resonances up to 2 GeV. One of the pur- 
poses of this study is to provide the results for the simulation codes of subthreshold kaon 
production in heavy ion collisions. This is the first, consistent study of the elementary 
kaon production reactions including both ttB and BB {B = N, A) collisions on the same 
footing. Comparisons are made between the calculated results and the existing semi- 
empirical parametrizations which are widely used for the simulation codes, as well as the 
experimental data. 



PACS: 24.10.Jv; 25.40.-h; 25.60.Dz; 25.70.-z; 25.70.Ef; 25.75.Dw 
Keywords: Kaon production; Elementary cross sections; Heavy ion collisions; Baryon 
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Strangeness production in heavy ion collisions is presently an issue of intense study. The 
enhancement of strangeness production in heavy ion collisions might indicate evidence for a new 
form of nuclear matter, such as a quark gluon plasma [0], 0- Particularly, kaon production in 
heavy ion collisions is considered a promising method not only of obtaining information about 
hot, dense nuclear matter, but also of determining the nuclear equation of state 0, [|, |j. 

For the study of kaon production in heavy ion collisions, the elementary cross sections in 
pion-nucleon and baryon-baryon collisions are some of the most basic and important ingredients 
for the microscopic calculation of the total kaon yield. However, the available experimental data 
for the baryon-baryon collision channels are very scarce (most exist for pp collision channels) 
and thus could not be parameterized in a satisfactory way for practical use. Indeed, there are 
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drastic differences between the commonly used, semi-empirical parametrizations for (the energy 
dependence of) the total kaon production cross sections in baryon-baryon collisions proposed 
by Randrup and Ko ||, and Schurmann and Zwermann [ff|. Thus, one has to parametrize these 
cross sections by relying on the theoretical investigations, which also describe the pion-baryon 
collision channels at the same time. 

Although simulation results for kaon production in heavy ion collisions suggest that the NA 
and A A collision channels give the dominant contribution for the total kaon yield || || |lO|j , 
the present situation regarding the theoretical investigations seems to be unsatisfactory. In the 
past, the reactions pp —>■ NYK (Y — A, E hyperons) were investigated within a One Boson 
Exchange Model (OBEM) by Ferrari Q, Yao [0, Wu and Ko [0, Laget and Deloff 



These calculations show quite contradictory contributions from the pion and kaon exchange 
processes, which can be ascribed to the different values of the coupling constants and cut-off 
parameters applied. Nevertheless, all of these calculations reproduce almost perfectly the small 
number of available experimental data points. This may imply the necessity of the constraints 
to fix the pion and kaon contributions consistently by investigating some other appropriate 
reactions. 

Recently, Li and Ko |16[, and Sibirtsev [T7| performed calculations of kaon production in 



heavy ion collisions and proton-nucleus collisions, respectively by using the elementary kaon 
production cross sections calculated using the OBEM. They adopted the results of the resonance 



model [|1| [19], g(J for the elementary cross sections nB — > YK in order to extract the B1B2 — > 
B3YK cross sections. There, an isospin average procedure is introduced for the kB — > YK 
cross sections, and the exchanged pion and kaon are treated on their mass shells. Furthermore, 



although the kaon exchange contribution was included in addition to the pion exchange |L6|, [L7 



the r] and p-meson exchange contributions were not studied. To be consistent with the resonance 
model, rj and p-meson exchanges should be included, since those resonances included in the 
model are observed to decay to the rjN and pN as well as the ttN (ttA) channels JH]]. Actually, 
as will be discussed later, the p-meson exchange gives a dominant contribution in the present 
calculations. 

In this article, we present the results for the energy dependence of the total cross sections 



pp — > NYK calculated in the resonance model [18, f9, Effl. For the other reactions, B\B 



B3YK (Bi, B2, -B3 = A, A), an extensive summary of the results of the study will be reported 



elsewhere [22 



The advantages of the present investigations are as follows. The resonance model [18, B5[ [20 



has already been applied to the study of the ttB — > YK (B = N, A) reactions. Although the 
model treatment is within the empirical tree level (observed widths for the resonances are 
used), it is fully relativistic, and includes all relevant baryon resonances up to 2 GeV, which are 
observed to decay to the hyperon and kaon channels. It turned out that the model was successful 



in explaining the abundance of available experimental data for the ttN — > YK reactions [24f|. 
This means that the richness of the data can provide a sound basis to the model because the 
parameters in the model are determined so as to reproduce the whole available data optimally. 
Thus, many of the effective interactions which are necessary to describe the B1B2 — > B^YK 
reactions are fixed within our model. Furthermore, we believe the present calculations are 
the first systematic calculations for both the 71 B — > YK and B1B2 — > B^YK (pp — > NYK) 
reactions. 

The processes relevant for the present calculations are depicted in Fig. [I]. Here B* stands 
for the baryon resonances which decay to the kaon-hyperon channels with their masses up to 
around 2 GeV @, namely A(1650)(f~), iV(1710)(i + ), iV(1720)(| + ) and A(1920)(| + ). In 



our model, kaons are always assumed to arise from these resonances at the same time as the 
hyperons. Once the effective Lagrangians and form factors for the vertices are fixed, most of 
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the relevant coupling constants can be calculated by using the experimental branching ratios. 
More detailed explanations of the model can be found in refs. |1| 0, pQf| . 

Effective Lagrangian densities necessary to evaluate the Feynman diagrams shown in Fig. [I 



are given by, 

C~kNN = -ignNNNj 5 TN ■ 7T, (1) 

Ari\W(1650) = -toV(1650) (N (1650)tN ■ 7T + NrN (1650) ■ 7T ) , (2) 

£*NN(17W) = -^JViV(mO)(iV(1710)75TiV-7f + iVf7 5 ^(l710)-7f), (3) 

C«nn { 1720) = ^ (1720) (N^1720)rN ■ d,n + NtN^1720) ■ d,n ) , (4) 

£^A(i920) = ^ (1920) (A"(1920)^V • + Nl V(1920) • ^tt) , (5) 

= -ig V NN^ 5 Nrj, (6) 

/Vivcmo) = -^iVJV(ino) (iV(1710) 75 A^ + iV7 5 ;V(1710)r/ ) , (7) 

/Wmo) = ^ 172 °) (^(1720)^ + ^(1720)^) , (8) 

A>w = ~QpNN i^N ^tN • Py, + iV ct^tN • d^p^J , (9) 

£ pAW( mo) = -^jv(mo)(jV(1710) 7 ^JV.^ + iVfy 4 iV(1710).^), (10) 

C pNN{172 o) = -ig P NN(i720) (^(1720)75^-^ + ^75^(1720)-^), (11) 

£XAAT(1650) = -^AAT(1650)(^(1650)AK + ^AiV(1650)), (12) 

ZW(mo) = -itoiV(mo) (iV(1710) 75 AK + ifA 75 iV(1710)) , (13) 

C KAN{im = 9KAN{1720) (ATM(1720)A^+ (d»K)AN»(1720)) , (14) 

^_ft"SAf(1710) = — ^^2^(1710) ^-^ r (1710)75T • I^ + KII -T7 5 iV(1710)) , (15) 

£jHW(irao) = ^ SjV(1720) f^(1720)f ■ £^# + (^)S • tW"(1720)) , (16) 

£kka(1920) = ^f 20) fc(1920)? ■ lud^K + (fyif) S ■ ?V(1920)) . (17) 

The notation and definitions appearing in the above equations are as follows: I is the transition 
operator defined by 

f M ,= £ (lil^M)el (18) 

£=±1,0 Z Z 



with M and /z being the third components of the isospin states, and r the Pauli matrices. 
N , N (1710) , N (1720) and A(1920) stand for the fields of the nucleon and baryon resonances. 
They are expressed by N = (p,n), similarly for the nucleon resonances, and A(1920) = 
(A(1920) ++ , A(1920) + , A(1920)°, A(1920)") in isospin space. The physical representations of 
the fields are, K T = (K+,K°) , K = (K-,^ ) , rr ± = (tti + m 2 )/V2, tt° = 7r 3 , similarly for 

the p-meson fields, and = (Si + i£ 2 )/v / 2, S° = S 3 , respectively, where the superscript T 
means the transposition operation. The meson fields are defined as annihilating (creating) the 
physical particle (anti-particle) states. For the propagators iSp(p) of the spin 1/2 and iG flu (p) 
of the spin 3/2 resonances we use: 

. „ , , 7 • V + m 



with 



(21) 



where m and r^ n// stand for the mass and full decay width of the corresponding resonances. 
For the form factors F(q) attached to the meson-resonance vertices, we use 

m = ^ (22) 

where g and A are the three momentum of the meson and cut-off parameter, respectively. The 
form factors, coupling constants and cut-off parameters are adopted from refs. ]TJ|, |HJ pOflf] 
For the cut-off parameters in the r)NB* and pNB* vertex form factors in pp (baryon-baryon) 
collisions which have not appeared in earlier work, we use the same values as those for the ttNB* 
vertices in order to reduce the number of new parameters. The other coupling constants, 
cut-off parameters and form factors are, as far as possible, taken from the Bonn nucleon- 
nucleon potential model [^3] (Model I in TABLE B.l). We use a dipole form factor and the 
tensor coupling constant is given by the ratio k = fpN n / QpN n = 6.1 for the pNN vertex. In 
order to show the structure of total amplitude in terms of the meson exchange and resonance 
intermediate states, we explicitly write down the contributions for the pp —> pAK + reaction as 
an example: 

M(pp-^pAK + ) = M(n, N (1Q50)) +M(n,N (1710)) +M(n,N (1720)) 

+ M(t], N(1710)) +M(v, N(1720)) + M(p, iV(1710)) + M(p, iV(1720)) 
+ exchange. (23) 



On the right hand side of eq. fl23|), the exchanged mesons and intermediate state resonances 
are written inside the brackets (c.f. Fig. [I]). Since there is no way to fix adequately the relative 
signs among the amplitudes, we simply neglect all the interference terms (including the exchange 
terms). For example, the total amplitude for pp —>■ pAK + given by eq. ( |2"3"D contains fourteen 
different amplitudes including the exchange amplitudes (diagram b) in Fig. |l|), which have 
relative minus signs to the corresponding direct amplitudes (diagram a) in Fig. |l|). Thus there 
arise 2 6 different possible relative sign combinations among them. Concerning the ineterference 
effects between the direct and exchange amplitudes, we comment as follows. At energy 1 
GeV above the threshold, we calculated the absolute ratios between the direct and exchange 
amplitudes for the p-meson exchange (which gives dominant contribution in our model) in order 
to make a rough estimate of the interefernce effects for the pp — * pAK + reaction. It turned 
out that the direct (exchange) process gave a dominant contribution when the final proton was 
scattered in the forward (backward) direction. The absolute ratios of the direct and exchange 
amplitudes were typically a few percent, when the absolute values of the amplitudes became 
large and gave dominant contributions to the cross sections. Furthermore, when the final proton 
stopped, a sum of the absolute values of the direct and the exchange amplitudes amounted to at 
most about a 10 percent of the forward (backward) scattering case. Although the interference 
effects between the direct and exchange processes could be important near the threshold, the 
effects seems to be small around 1 GeV above the threshold, where the new parameters of the 

lr The coupling constants for 5^-£jv(mo) anc ^ Sif£JV(i720) m ^ ne tables of refs. JlJ, |l^] were listed wrongly, 
although the numerical calculations were performed by using the correct values as given in this article. 
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Table 1: Coupling constants and cut-off parameters, k = fpNN 1 9 p nn = 6.1 for the pNN 
tensor coupling is used. The confidence levels of the resonances listed below are, N(1650)****, 
N(1710)***, N(1720)**** and A(1920)*** 
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model were fitted to the data. The effects of interference terms were studied in refs. [ 18| for 
the collision channels, and it turned out that their effects to the shape of the differential 
cross sections were appreciable, but their effects to the total cross sections were small. 

For completeness, we give the relations between the K + and K° production cross sections 
treated in this article as follows: 

a (pp^pAK + ) = a(nn — > nAK°), (24) 

a(nn -> n5TAT + ) = a(pp -> p£ + AT ), (25) 

a(pp ^pY?K + ) = a(nn^nE°K°), (26) 

a(pp -> nZ + K + ) = a(nn^nS~K°). (27) 

Here it should be mentioned that the p meson exchange turned out to give a dominant 
contribution in the present calculations, as mentioned before. In order to reproduce all the 
available data in a satisfactory manner, the coupling constant g P NN and cut-off parameter A p 
were varied. The dependence on these quantities is discussed below. 

In Fig. 0, the dependence on the coupling constant g P NN of the total cross section pp — > 
pAK + is shown. Here, s 1//2 is the invariant collision energy in the proton-proton center-of- 

1 /2 

momentum system, and s = + ray + is the threshold energy, with m^, my an d 
rriK being the masses of the nucleon, hyperon and kaon, respectively. After fixing the cut-off 
parameter A p to a specific value, the dependence on the coupling constant is small. 

Sensitivity of the total cross section pp — ► pAK + to the cut-off parameter values A p are 
shown in Fig. |3] with the experimental data |24j]. The dependence on the cut-off parameter A p 
is rather large in the present calculation, after fixing the coupling constant q p nn (and k=6.1 
for the tensor coupling constant) to adequate, non-small values. 

After all, the best values obtained to reproduce the experimental data are, g 2 pNN /A'K = 0.74 
and A p = 920 MeV. We summarize in Table [l] all the relevant coupling constants and cut-off 
parameters determined and used for the calculations in our model. 

In Fig. |] we show the calculated energy dependence of the total cross sections pp — > pY, + K°, 
pp — > pTPK + and pp — > nH + K + in comparison with the experimental data pi. The variations 
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of the experimental data are rather large, but our model can simultaneously reproduce these 
data well. Remember that we do not have so much freedom in the number of parameters to be 
varied. 

Now we compare our results with some of the existing semi-empirical parametrizations 
for the energy dependence of the total cross section pp — > pAK + . In Fig. [5], the dashed- 
dotted, the dotted, the dashed and the solid lines stand for the parametrizations of Randrup 
and Ko ||, Schiirmann and Zwermann 0, OBEM results of Sibirtsev [17], and our results, 



respectively. They illustrate quite different energy dependence and strongly disagree when the 
collision energy, s 1//2 , is close to the reaction threshold, except that the results of Sibirtsev [17 



and the present calculations show similar behaviour. We should mention here that Li and Ko 



Tq| obtained similar results to ours by applying kaon and pion exchanges. There, the pion and 
kaon exchanges give almost the same contributions for the pp — > pAK + total cross sections, 
while the pion exchange alone gives almost sufficient contribution for the pp — > pT*K cross 
sections. In our model, however, kaon exchange does not enter the calculation, since kaons are 
assumed to appear together with hyperons through the intermediate state baryon resonances. 
Furthermore, the exchanged mesons in our model, 7T, rj and p mesons are adopted based on the 
observed decay channels of these baryon resonances. However, the contribution of the rj meson 
exchange in our model is small. 

We believe the strangeness production in pp collisions, which is currently being performed at 
COSY-Julich [p5|| , will provide important data to improve our understanding of the elementary 
kaon production mechanism. 

In summary, we have presented new calculations for the energy dependence of the total 
kaon production cross sections in proton-proton collisions using the resonance model. The re- 
sults reported in this article are part of the systematic studies for both the pion-baryon and 
baryon-baryon collision channels. The investigations for free space will be complete in the near 
future P2[ . The next task to be done is to incorporate the effects of the medium on these ele- 
mentary kaon production cross sections, which can be implimented in a unified manner based 
on the same model. Then, we hope that those in-medium modified cross sections can be applied 
to more realistic investigations of kaon production in heavy ion and proton-nucleus collisions. 
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Figure 1: Processes contribute to kaon production in pp collisions. In figure, B* stands for the 
resonances which are observed to decay to kaon (K) and hyperons (Y = A, E) up to their masses 
around 2 GeV. They are, 7V(1650)(§~), iV(I7I0)(± + ), iV(I720)(§ + ) and A(I920)(§ + ). The 
exchanged mesons 7r, i] and p are experimentally observed in these resonance decay channels. 
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Figure 2: Dependence on the coupling constant g P NN of the total cross section pp —>■ pAK + . The 
dotted, the solid, and the dashed lines show the results obtained with the coupling constant 
0.84, ff^ VAr /47r=0.74, and gp NN /4n=0.70 with the cut-off parameter A p = 



values g 2 pNN /4:7i 



920 MeV, respectively. The dots show experimental data p3 with error bars. 




pAK + to the cut-off parameters A p with 
the dashed and the solid lines are the results obtained with 
the cut-off parameter values A p = 1 000 MeV, A p = 940 MeV and A p = 920 MeV, respectively. 
The coupling constant value is fixed to g 2 NN /An 



Figure 3: Sensitivities of the total cross section pp 
experimental data [24]. The dotted 



0.74. 
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Figure 4: Energy dependence of the total cross sections for the pp — > pTi + K°, pp — > pYPK + 
and pp — > nYi + K + reactions. The dots stand for the experimental data with error bars. 
The solid lines show our results. 




Figure 5: Comparisons between the existing semi-empirical parametrizations and our results for 
the energy dependence of the total cross section pp — > pAK + . The dotted, the dashed-dotted, 
the dashed, and the solid lines stand for the parametrizations of Schurmann and Zwermann |7J , 
Randrup and Ko [0], OBEM results of Sibirtsev fl7j, and our results, respectively. 
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